A Comprehensive Allele Specific Expression Resource for the Equine Transcriptome.

Background: Allele-specific expression (ASE) analysis provides a nuanced view of cis-regulatory mechanisms affecting gene expression. Results: An equine ASE analysis was performed, using integrated Iso-seq and short-read RNA sequencing data from four healthy Thoroughbreds (2 mares and 2 stallions) across 9 tissues from the Functional Annotation of Animal Genomes (FAANG) project. Allele expression was quantified by haplotypes from long-read data, with 42,900 allele expression events compared. Within these events, 635 (1.48%) demonstrated ASE, with liver tissue containing the highest proportion. Genetic variants within ASE events were in histone modified regions 64.2% of the time. Validation of allele-specific variants, using a set of 66 equine liver samples from multiple breeds, confirmed that 97% of variants demonstrated ASE. Conclusions: This valuable publicly accessible resource is poised to facilitate investigations into regulatory variation in equine tissues. Our results highlight the tissue-specific nature of allelic imbalance in the equine genome.

A Comprehensive Allele Specific Expression Resource for the Equine Transcriptome.

Background: Allele-specific expression (ASE) analysis provides a nuanced view of cis-regulatory mechanisms affecting gene expression. Results: In this work, we introduce and highlight the significance of an equine ASE analysis, containing integrated long- and short-read RNA sequencing data, along with insight from histone modification data, from four healthy Thoroughbreds (2 mares and 2 stallions) across 9 tissues. Conclusions: This valuable publicly accessible resource is poised to facilitate investigations into regulatory variation in equine tissues and foster a deeper understanding of the impact of allelic imbalance in equine health and disease at the molecular level.

Genetic polymorphisms in vitamin E transport genes as determinants for risk of equine neuroaxonal dystrophy.

BACKGROUND: Equine neuroaxonal dystrophy/equine degenerative myeloencephalopathy (eNAD/EDM) is an inherited neurodegenerative disorder associated with vitamin E deficiency. In humans, polymorphisms in genes involved in vitamin E uptake and distribution determines individual vitamin E requirements. HYPOTHESIS/OBJECTIVES: Genetic polymorphisms in genes involved in vitamin E metabolism would be associated with an increased risk of eNAD/EDM in Quarter Horses (QHs). ANIMALS: Whole-genome sequencing: eNAD/EDM affected (n = 9, postmortem [PM]-confirmed) and control (n = 32) QHs. VALIDATION: eNAD/EDM affected (n = 39, 23-PM confirmed) and control (n = 68, 7-PM confirmed) QHs. Allele frequency (AF): Publicly available data from 504 horses across 47 breeds. METHODS: Retrospective, case control study. Whole-genome sequencing was performed and genetic variants identified within 28 vitamin E candidate genes. These variants were subsequently genotyped in the validation cohort. RESULTS: Thirty-nine confirmed variants in 15 vitamin E candidate genes were significantly associated with eNAD/EDM (P < .01). In the validation cohort, 2 intronic CD36 variants (chr4:726485 and chr4:731082) were significantly associated with eNAD/EDM in clinical (P = 2.78 × 10-4 and P = 4 × 10-4 , respectively) and PM-confirmed cases (P = 6.32 × 10-6 and 1.04 × 10-5 , respectively). Despite the significant association, variant AFs were low in the postmortem-confirmed eNAD/EDM cases (0.22-0.26). In publicly available equine genomes, AFs ranged from 0.06 to 0.1. CONCLUSIONS AND CLINICAL IMPORTANCE: Many PM-confirmed cases of eNAD/EDM were wild-type for the 2 intronic CD36 SNPs, suggesting either a false positive association or genetic heterogeneity of eNAD/EDM within the QH breed.

The localization of centromere protein A is conserved among tissues.

Centromeres are epigenetically specified by the histone H3 variant CENP-A. Although mammalian centromeres are typically associated with satellite DNA, we previously demonstrated that the centromere of horse chromosome 11 (ECA11) is completely devoid of satellite DNA. We also showed that the localization of its CENP-A binding domain is not fixed but slides within an about 500 kb region in different individuals, giving rise to positional alleles. These epialleles are inherited as Mendelian traits but their position can move in one generation. It is still unknown whether centromere sliding occurs during meiosis or during development. Here, we first improve the sequence of the ECA11 centromeric region in the EquCab3.0 assembly. Then, to test whether centromere sliding may occur during development, we map the CENP-A binding domains of ECA11 using ChIP-seq in five tissues of different embryonic origin from the four horses of the equine FAANG (Functional Annotation of ANimal Genomes) consortium. Our results demonstrate that the centromere is localized in the same region in all tissues, suggesting that the position of the centromeric domain is maintained during development.

Functional annotation of the animal genomes: An integrated annotation resource for the horse.

The genomic sequence of the horse has been available since 2009, providing critical resources for discovering important genomic variants regarding both animal health and population structures. However, to fully understand the functional implications of these variants, detailed annotation of the horse genome is required. Due to the limited availability of functional data for the equine genome, as well as the technical limitations of short-read RNA-seq, existing annotation of the equine genome contains limited information about important aspects of gene regulation, such as alternate isoforms and regulatory elements, which are either not transcribed or transcribed at a very low level. To solve above problems, the Functional Annotation of the Animal Genomes (FAANG) project proposed a systemic approach to tissue collection, phenotyping, and data generation, adopting the blueprint laid out by the Encyclopedia of DNA Elements (ENCODE) project. Here we detail the first comprehensive overview of gene expression and regulation in the horse, presenting 39,625 novel transcripts, 84,613 candidate cis-regulatory elements (CRE) and their target genes, 332,115 open chromatin regions genome wide across a diverse set of tissues. We showed substantial concordance between chromatin accessibility, chromatin states in different genic features and gene expression. This comprehensive and expanded set of genomics resources will provide the equine research community ample opportunities for studies of complex traits in the horse.

Investigation of high gamma-glutamyltransferase syndrome in California Thoroughbred racehorses.

BACKGROUND: Increases in serum gamma-glutamyltransferase (GGT) activity have been reported in Thoroughbred (TB) racehorses and associated with maladaptation to training but the underlying etiology remains unknown. HYPOTHESIS/OBJECTIVES: Classify the etiology of high GGT syndrome in racing TBs by assessment of pancreatic enzymes, vitamin E concentrations, and both a candidate gene and whole genome association study. We hypothesized that a genetic variant resulting in antioxidant insufficiency or pancreatic dysfunction would be responsible for high GGT syndrome in TBs. ANIMALS: A total of 138 California racing TBs. Amylase: n = 31 affected (serum GGT activity ≥60 IU/L), n = 52 control (serum GGT activity <40 IU/L). Lipase: n = 19 affected, n = 35 control. Serum α-tocopherol concentrations: n = 32 affected, n = 46 control. Genome-wide association study (GWAS): 36 affected, 58 control. Whole genome sequencing: n = 5 affected, n = 5 control. METHODS: Biochemical and vitamin analytes were compared among cohorts. A GWAS was performed and a subset of TBs underwent whole genome sequencing to interrogate candidate genes and positional genetic regions. RESULTS: Serum lipase and amylase activity and α-tocopherol concentrations did not differ between groups. No genetic variants were identified in 2 candidate genes (UGT1A1 and GGT1) that associated with the phenotype. Four single nucleotide polymorphisms (SNPs) approached a suggestive association with the phenotype (P = 2.15 × 10-5 ), defining a 100 kb region on chromosome 5 surrounding cluster of differentiation 1a (CD1A1), a transmembrane gene related to the major histocompatibility complex. CONCLUSIONS AND CLINICAL IMPORTANCE: An underlying genetic etiology may exist for high GGT syndrome in racing TBs, similar to genetic disorders in humans.

DNA methylation aging and transcriptomic studies in horses.

Cytosine methylation patterns have not yet been thoroughly studied in horses. Here, we profile n = 333 samples from 42 horse tissue types at loci that are highly conserved between mammalian species using a custom array (HorvathMammalMethylChip40). Using the blood and liver tissues from horses, we develop five epigenetic aging clocks: a multi-tissue clock, a blood clock, a liver clock and two dual-species clocks that apply to both horses and humans. In addition, using blood methylation data from three additional equid species (plains zebra, Grevy's zebras and Somali asses), we develop another clock that applies across all equid species. Castration does not significantly impact the epigenetic aging rate of blood or liver samples from horses. Methylation and RNA data from the same tissues define the relationship between methylation and RNA expression across horse tissues. We expect that the multi-tissue atlas will become a valuable resource.

Decoding the Equine Genome: Lessons from ENCODE.

The horse reference genome assemblies, EquCab2.0 and EquCab3.0, have enabled great advancements in the equine genomics field, from tools to novel discoveries. However, significant gaps of knowledge regarding genome function remain, hindering the study of complex traits in horses. In an effort to address these gaps and with inspiration from the Encyclopedia of DNA Elements (ENCODE) project, the equine Functional Annotation of Animal Genome (FAANG) initiative was proposed to bridge the gap between genome and gene expression, providing further insights into functional regulation within the horse genome. Three years after launching the initiative, the equine FAANG group has generated data from more than 400 experiments using over 50 tissues, targeting a variety of regulatory features of the equine genome. In this review, we examine how valuable lessons learned from the ENCODE project informed our decisions in the equine FAANG project. We report the current state of the equine FAANG project and discuss how FAANG can serve as a template for future expansion of functional annotation in the equine genome and be used as a reference for studies of complex traits in horse. A well-annotated reference functional atlas will also help advance equine genetics in the pan-genome and precision medicine era.

Successful ATAC-Seq From Snap-Frozen Equine Tissues.

An assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) has become an increasingly popular method to assess genome-wide chromatin accessibility in isolated nuclei from fresh tissues. However, many biobanks contain only snap-frozen tissue samples. While ATAC-seq has been applied to frozen brain tissues in human, its applicability in a wide variety of tissues in horse remains unclear. The Functional Annotation of Animal Genome (FAANG) project is an international collaboration aimed to provide high quality functional annotation of animal genomes. The equine FAANG initiative has generated a biobank of over 80 tissues from two reference female animals and experiments to begin to characterize tissue specificity of genome function for prioritized tissues have been performed. Due to the logistics of tissue collection and storage, extracting nuclei from a large number of tissues for ATAC-seq at the time of collection is not always practical. To assess the feasibility of using stored frozen tissues for ATAC-seq and to provide a guideline for the equine FAANG project, we compared ATAC-seq results from nuclei isolated from frozen tissue to cryopreserved nuclei (CN) isolated at the time of tissue harvest in liver, a highly cellular homogenous tissue, and lamina, a relatively acellular tissue unique to the horse. We identified 20,000-33,000 accessible chromatin regions in lamina and 22-61,000 in liver, with consistently more peaks identified using CN isolated at time of tissue collection. Our results suggest that frozen tissues are an acceptable substitute when CN are not available. For more challenging tissues such as lamina, nuclei extraction at the time of tissue collection is still preferred for optimal results. Therefore, tissue type and accessibility to intact nuclei should be considered when designing ATAC-seq experiments.

Increased α-tocopherol metabolism in horses with equine neuroaxonal dystrophy.

BACKGROUND: Equine neuroaxonal dystrophy/equine degenerative myeloencephalopathy (eNAD/EDM) is an inherited neurodegenerative disorder associated with a vitamin E deficiency within the first year of life. Vitamin E consists of 8 isoforms metabolized by the CYP4F2 enzyme. No antemortem diagnostic test currently exists for eNAD/EDM. HYPOTHESIS/OBJECTIVES: Based on the association of α-tocopherol deficiency with the development of eNAD/EDM, we hypothesized that the rate of α-tocopherol, but not γ-tocopherol or tocotrienol metabolism, would be increased in eNAD/EDM-affected horses. ANIMALS: Vitamin E metabolism: Proof of concept (POC) study; eNAD/EDM-affected (n = 5) and control (n = 6) horses. Validation study: eNAD/EDM-affected Quarter Horses (QHs; n = 6), cervical vertebral compressive myelopathy affected (n = 6) horses and control (n = 29) horses. CYP4F2 expression and copy number: eNAD/EDM-affected (n = 12) and age- and sex-matched control (n = 11-12) horses. METHODS: The rates of α-tocopherol/tocotrienol and γ-tocopherol/tocotrienol metabolism were assessed in equine serum (POC and validation) and urine (POC only) using liquid chromatography tandem mass spectrometry (LC-MS/MS). Quantitative reverse-transcriptase PCR (qRT-PCR) and droplet digital (dd)-PCR were used to assay expression and genomic copy number of a CYP4F2 equine ortholog. RESULTS: Metabolic rate of α-tocopherol was increased in eNAD/EDM horses (POC,P < .0001; validation, P = .03), with no difference in the metabolic rate of γ-tocopherol. Horses with eNAD/EDM had increased expression of the CYP4F2 equine orthologue (P = .02) but no differences in copy number. CONCLUSIONS AND CLINICAL IMPORTANCE: Increased α-tocopherol metabolism in eNAD/EDM-affected QHs provides novel insight into alterations in vitamin E processing in eNAD/EDM and highlights the need for high-dose supplementation to prevent the clinical phenotype in genetically susceptible horses.

Generation of a Biobank From Two Adult Thoroughbred Stallions for the Functional Annotation of Animal Genomes Initiative.

Following the successful creation of a biobank from two adult Thoroughbred mares, this study aimed to recapitulate sample collection in two adult Thoroughbred stallions as part of the Functional Annotation of the Animal Genome (FAANG) initiative. Both stallions underwent thorough physical, lameness, neurologic, and ophthalmic (including electroretinography) examinations prior to humane euthanasia. Epididymal sperm was recovered from both stallions immediately postmortem and cryopreserved. Aseptically collected full thickness skin biopsies were used to isolate, culture and cryopreserve dermal fibroblasts. Serum, plasma, cerebrospinal fluid, urine, and gastrointestinal content from various locations were collected and cryopreserved. Under guidance of a board-certified veterinary anatomic pathologist, 102 representative tissue samples were collected from both horses. Whole tissue samples were flash-frozen and prioritized tissues had nuclei isolated and cryopreserved. Spatially contemporaneous samples of each tissue were submitted for histologic examination. Antemortem and gross pathologic examination revealed mild abnormalities in both stallions. One stallion (ECA_UCD_AH3) had unilateral thoracic limb lameness and bilateral chorioretinal scars. The second stallion (ECA_UCD_AH4) had subtle symmetrical pelvic limb ataxia, symmetrical prostatomegally, and moderate gastrointestinal nematodiasis. DNA from each was whole-genome sequenced and genotyped using the GGP Equine 70K SNP array. The genomic resources and banked biological samples from these animals augments the existing resource available to the equine genomics community. Importantly we may now improve the resolution of tissue-specific gene regulation as affected by sex, as well as add sex-specific tissues and gametes.

Commercial genetic testing for type 2 polysaccharide storage myopathy and myofibrillar myopathy does not correspond to a histopathological diagnosis.

BACKGROUND: Commercial genetic tests for type 2 polysaccharide storage myopathy (PSSM2) and myofibrillar myopathy (MFM) have not been validated by peer-review, and formal regulation of veterinary genetic testing is lacking. OBJECTIVES: To compare genotype and allele frequencies of commercial test variants (P variants) in MYOT (P2; rs1138656462), FLNC (P3a; rs1139799323), FLNC (P3b; rs1142918816) and MYOZ3 (P4; rs1142544043) between Warmblood (WB) and Arabian (AR) horses diagnosed with PSSM2/MFM by muscle histopathology, and phenotyped breed-matched controls. To quantify variant frequency in public repositories of ancient and modern horse breeds. STUDY DESIGN: Cross sectional using archived clinical material and publicly available data. METHODS: We studied 54 control-WB, 68 PSSM2/MFM-WB, 30 control-AR, 30 PSSM2/MFM-AR and 205 public genotypes. Variants were genotyped by pyrosequencing archived DNA. Genotype and allele frequency, and number of variant alleles or loci were compared within breed between controls, PSSM2/MFM combined and MFM or PSSM2 horses considered separately using additive/genotypic and dominant models (Fisher's exact tests). Variant frequencies in modern, early domestic and Przewalski horses were determined from a public data repository. RESULTS: There was no significant association between any P locus and a histopathological diagnosis of PSSM2/MFM, and no difference between control and myopathic horses in total loci with alternative alleles, or total alternate alleles when PSSM2/MFM was considered combined or separately as PSSM2 or MFM. For all tests, sensitivity was <0.33. Allele frequencies in WB (controls/cases) were: 8%/15% (P2), 5%/6% (P3a/b) and 9%/13% (P4); in AR, frequencies were: 12%/17% (P2), 2%/2% (P3a/b) and 7%/12% (P4). All P variants were present in early domestic (400- to 5500-year-old) horses and P2 present in the Przewalski. CONCLUSIONS: Because of the lack of significant association between a histopathological diagnosis of PSSM2 or MFM and the commercial genetic test variants P2, P3 and P4 in WB and AR, we cannot recommend the use of these variant genotypes for selection and breeding, prepurchase examination or diagnosis of a myopathy.

Mandibulofacial Dysostosis Attributed to a Recessive Mutation of CYP26C1 in Hereford Cattle.

In spring 2020, six Hereford calves presented with congenital facial deformities attributed to a condition we termed mandibulofacial dysostosis (MD). Affected calves shared hallmark features of a variably shortened and/or asymmetric lower mandible and bilateral skin tags present 2-10 cm caudal to the commissure of the lips. Pedigree analysis revealed a single common ancestor shared by the sire and dam of each affected calf. Whole-genome sequencing (WGS) of 20 animals led to the discovery of a variant (Chr26 g. 14404993T>C) in Exon 3 of CYP26C1 associated with MD. This missense mutation (p.L188P), is located in an α helix of the protein, which the identified amino acid substitution is predicted to break. The implication of this mutation was further validated through genotyping 2 additional affected calves, 760 other Herefords, and by evaluation of available WGS data from over 2500 other individuals. Only the affected individuals were homozygous for the variant and all heterozygotes had at least one pedigree tie to the suspect founder. CYP26C1 plays a vital role in tissue-specific regulation of retinoic acid (RA) during embryonic development. Dysregulation of RA can result in teratogenesis by altering the endothelin-1 signaling pathway affecting the expression of Dlx genes, critical to mandibulofacial development. We postulate that this recessive missense mutation in CYP26C1 impacts the catalytic activity of the encoded enzyme, leading to excess RA resulting in the observed MD phenotype.

A nonsense variant in Rap Guanine Nucleotide Exchange Factor 5 (RAPGEF5) is associated with equine familial isolated hypoparathyroidism in Thoroughbred foals.

Idiopathic hypocalcemia in Thoroughbred (TB) foals causes tetany and seizures and is invariably fatal. Based upon the similarity of this disease with human familial hypoparathyroidism and occurrence only in the TB breed, we conducted a genetic investigation on two affected TB foals. Familial hypoparathyroidism was identified, and pedigree analysis suggested an autosomal recessive (AR) mode of inheritance. We performed whole-genome sequencing of the two foals, their unaffected dams and four unaffected, unrelated TB horses. Both homozygosity mapping and an association analysis were used to prioritize potential genetic variants. Of the 2,808 variants that significantly associated with the phenotype using an AR mode of inheritance (P<0.02) and located within a region of homozygosity, 1,507 (54%) were located in a 9.7 Mb region on chr4 (44.9-54.6 Mb). Within this region, a nonsense variant (RAPGEF5 c.2624C>A,p.Ser875*) was significantly associated with the hypoparathyroid phenotype (Pallelic = 0.008). Affected foals were homozygous for the variant, with two additional affected foals subsequently confirmed in 2019. Necropsies of all affected foals failed to identify any histologically normal parathyroid glands. Because the nonsense mutation in RAPGEF5 was near the C-terminal end of the protein, the impact on protein function was unclear. Therefore, we tested the variant in our Xenopus overexpression model and demonstrated RAPGEF5 loss-of-function. This RAPGEF5 variant represents the first genetic variant for hypoparathyroidism identified in any domestic animal species.

Genome-Wide Association Study and Subsequent Exclusion of ATCAY as a Candidate Gene Involved in Equine Neuroaxonal Dystrophy Using Two Animal Models.

Equine neuroaxonal dystrophy/equine degenerative myeloencephalopathy (eNAD/EDM) is an inherited neurodegenerative disorder of unknown etiology. Clinical signs of neurological deficits develop within the first year of life in vitamin E (vitE) deficient horses. A genome-wide association study (GWAS) was carried out using 670,000 SNP markers in 27 case and 42 control Quarter Horses. Two markers, encompassing a 2.5 Mb region on ECA7, were associated with the phenotype (p = 2.05 × 10-7 and 4.72 × 10-6). Within this region, caytaxin (ATCAY) was identified as a candidate gene due to its known role in Cayman Ataxia and ataxic/dystonic phenotypes in mouse models. Whole-genome sequence data in four eNAD/EDM and five unaffected horses identified 199 associated variants within the ECA7 region. MassARRAY® genotyping was performed on these variants within the GWAS population. The three variants within ATCAY were not concordant with the disease phenotype. No difference in expression or alternative splicing was identified using qRT-PCR in brainstem across the ATCAY transcript. Atcayji-hes mice were then used to conduct functional analysis in a second animal model. Histologic lesions were not identified in the central nervous system of Atcayji-hes mice. Additionally, supplementation of homozygous Atcayji-hes mice with 600 IU/day of dl-α-tocopheryl acetate (vitE) during gestation, lactation, and adulthood did not improve the phenotype. ATCAY has therefore been excluded as a candidate gene for eNAD/EDM.

A missense mutation in MYH1 is associated with susceptibility to immune-mediated myositis in Quarter Horses.

BACKGROUND: The cause of immune-mediated myositis (IMM), characterized by recurrent, rapid-onset muscle atrophy in Quarter Horses (QH), is unknown. The histopathologic hallmark of IMM is lymphocytic infiltration of myofibers. The purpose of this study was to identify putative functional variants associated with equine IMM. METHODS: A genome-wide association (GWA) study was performed on 36 IMM QHs and 54 breed matched unaffected QHs from the same environment using the Equine SNP50 and SNP70 genotyping arrays. RESULTS: A mixed model analysis identified nine SNPs within a ~ 2.87 Mb region on chr11 that were significantly (Punadjusted < 1.4 × 10- 6) associated with the IMM phenotype. Associated haplotypes within this region encompassed 38 annotated genes, including four myosin genes (MYH1, MYH2, MYH3, and MYH13). Whole genome sequencing of four IMM and four unaffected QHs identified a single segregating nonsynonymous E321G mutation in MYH1 encoding myosin heavy chain 2X. Genotyping of additional 35 IMM and 22 unaffected QHs confirmed an association (P = 2.9 × 10- 5), and the putative mutation was absent in 175 horses from 21 non-QH breeds. Lymphocytic infiltrates occurred in type 2X myofibers and the proportion of 2X fibers was decreased in the presence of inflammation. Protein modeling and contact/stability analysis identified 14 residues affected by the mutation which significantly decreased stability. CONCLUSIONS: We conclude that a mutation in MYH1 is highly associated with susceptibility to the IMM phenotype in QH-related breeds. This is the first report of a mutation in MYH1 and the first link between a skeletal muscle myosin mutation and autoimmune disease.